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Abstract — In this paper we consider stationary Forced KdV 
equation with positive Forcing term. The supercritical solitary 
wave solutions of the stationary Forced KdV equation are 
obtained. In order to obtain the solutions the domain of the 
problem has been divided in to three parts; the left, the middle 
and the right parts. The solution on the left and the right parts 
are obtained by an analytical method. The solution on the 
middle part is expressed in the terms of Weierstarss elliptic 
function. We have designed computer programs using 
Mathematica to produce the solutions. The complete solution 
was found by matching the solutions of all the three parts. We 
have found out that there are four different solutions according 
to the values of the phase shift. Only one solution is positive. 
Further research can be carried out for negative forcing terms. 

Index Terms — Stationary Forced KdV Equation; 
Supercritical solution;, positive Forcin 


I. INTRODUCTION 

The forced KdV equation is a first -order approximation of a 
long nonlinear surface wave in a channel flow of an in viscid 
fluid of constant density over a bump or a dent. This equation 
is normally written in the form (Gong and Shen, 1994) 

r] t + Arj x + 2arjt] x + Prj m = f x [x) -ao<x<ao f >0 ( 1 ) 

where A > 0, (X < 0 and J3 <0 are constant and f(x) is a 
given function called the forcing term which is differentiable 
and has a compact support i.e. it is nonzero only in a closed 
bounded set. Equation (1) was first derived by Akylas (1984) 
and in an asymptotically reduced result from Euler equations 
of fluid motion and the corresponding boundary conditions. 
The unknown TJ[X,t ) represents the first order elevation of 
the free surface of the fluid from it is equilibrium level. The 
forcing function f(x) is due to the bottom topography of the 
fluid domain (such as a bump or dent) or due to an external 
pressure on the free surface such as the wind stress on the 
surface of an ocean. In the absence of the forcing function 
term i.e. f(x) = 0 equation (1) becomes the familiar Korteweg 
deVries(KdV) equation 

rjt + Xrix + 'iGtTjrix + [Srj^ = 0 -co<x<oo ?>0 ( 2 ) 

when T] t = 0 equation (1) becomes the equation of our 
concern, the stationary forced (sfKdV) equation, 
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Arj x + 2arjr] x + = f x {x) -oo<;c<ao (3) 

Equation (3) was first derived by Shen (1993). The solutions 
of equation (3) are categorized according to the value of 
A (Upstream flow velocity) as follows (Shen, 1993). 

(i) Supercritical stationary waves. They occur only when A is 
positive and sufficiently large. 

(ii) Subcritical stationary waves. They occur only when A is 
negative and sufficiently small. 

(hi) Unsteady periodic soliton radiation. This solution 
appears when |/l| is small. Such a solution is called the 
transcritical solution. 

n. Problem Statement 

In this paper we study supercritical solitary wave solutions of 
equation (3) with positive forcing terms in the rectangular 
bump. In this work we shall take the value of A equal to 3 
and the length of rectangular bump equals 2 in order to 
compare with the previous results. 

Solitary waves refer to any surface wave profile that dies out 
at infinity which means that the free surface elevation r] has 
the property 

77 -> 0, rj x -A- 0 and rj xi — >0 as x —> ±00 (4) 

Integrating equation (3) once with respect to the independent 
variable from -00 to x, we have 

Arj + arf + prj^ = /(*), - oo < x < 00 , A > 0 (5) 

with 77 ^ 0, and tj x —» 0 as x —> ±00 

In this paper we shall study equation (5) and (6) with the 
following conditions 

3 1 

(i) a = — <0, J3 = — and A = 3. 

4 H 6 

[l, ~-<x<- 

(ii) f{x) = 2 2 

[0 otherwise 

where a is positive constant representing the length of a 
rectangular bump. 
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In order to obtain the solution of equation (5) and (6) the 
domain of the problem has been divided in three parts. The 
left (x < -a/2), the middle (-a/2 < x < a/2) and the right parts (x 
> a/2). The complete solution was found by matching the 
solutions of all parts. Section 3 and Section 4 describes the 
method of finding solution on the left and the right part 
respectively. Section 5 describes the solution in the middle 
part. In Section 6 we determine the phase shift L a which 
depends on the forcing terms and the continuity conditions 
between each adjacent part. In Section 7 we combine the 
solutions obtained in Section 4 and Section 5. The 
conclusions of this paper are summarized in Section 8. 


III. Left-Side Solution 

In this Section we need to solve equations (5) and (6) in the 
region of zero forcing term by using analytical method. 
Meanwhile, the solution in the region of non zero forcing term 
will be discussed later in Section 5. 

Outside the region of the forcing term, equation (5) and (6) 
becomes 


Atj + arp +Prj xx = 0 ^<x<-^, A>Q (7) 

7^0, and rj x . 0 —> 0 as x —» ±oo (8) 


We need to find the solution in the negative region (x < -a/2) outside the rectangular bump. In this region 
equations (7) and (8) are written as 

Arj + arj 1 + /3ij" =0 — co<x<—^, A>0 (9) 

7 ^ 0 , and rj x —» 0 —» 0 as x —> -oo (10) 
where equation (10) is the half negative solitary wave conditions. We start by rearranging equation (9) in 
the form 


Ptj" -—Atj - ar] 1 -qo<x<-— , A> 0 (11) 


/ dr/ 

where 7 =- 

dx 


The procedure we adopt here in order to integrate equation (11) is that we write 7 as 

/ drj' 


a = 7 ~r~ (12) 


drj 


Substituting equation (12) into equation (11) we obtain 


pTj'drj 1 = {^Ar]-arf)flri 


Equation (13) is simply a first order separable equation which can be integrated easily to get 

P ( / )2 A 2 a 3 ^ 

ffo) = -2 17 y 17 +c 


(13) 

(14) 


Multiplying both sides of equation(14) by — we then obtain 

3 + 2 c 

2a 2 a a 


(15) 


where — is positive constant and C is the integration constant. 

a 

By the half negative solitary wave conditions (10), equation (15) has a real double root, r o and r o is 
smaller than the third real root r. Then equation (15) can be written as 

f(/r- {r-rih-r,\r.<r (16) 
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If we now let 



then equation (16) can be written a; 




s l =r-r o >0. (17) 


Equation (17) can be integrated to yield 


Hence, the actual solution on the left side is 


t* 3S 1 A 

v = s 1 seen J—£ 


" w= ~^ sech2 M (x ~ Lj 


where L a is the phase shift to be determined in Section 6. 
All graphs of these solutions are shown in Figure 1. In 
Section 6 we shall show how to calculate L a . For the purpose 
of illustration we now use the values of L a obtained in 
Section 6 to produce the appropriate solution of the sfKdV on 
the left 


The graphs shown in Figure 1 are the solution of the sfKdV 
equation on the left side when A = 3 and a = 2. The values of 
L a used for this purpose are -0.20752, -0.797065, -1.792481 
and -2.248243 for graph (a), (b), (c) and (d) respectively. 



Figure 1: Feft-Side Solution 

IV. The Right-Side Solution 

a 

In this section we determine the solution on the right side (X > — ) that out side the rectangular bump. In 
this region equations (7) and (8) become 
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Xr] + ar/ 2 + fir] 11 = 0, x <co, X > 0 (20) 

77 —> 0 , 77 ^ —> 0 as x— »+oo ( 21 ) 

where equation (21) in normally termed as the positive half solitary wave conditions. 

The different between the left and the right region is simply on the solitary wave conditions. The same 
procedure as used for the left side can be applied to find the solution in the right and we thus obtain 

flto 4 * (22) 

where the phase shift L x = —L 0 . All graphs of this solution are shown in Figure 2. 

The graphs shown in Figure 2 are the solution of the sfKdV equation on the right side when X = 3 and a = 
2. The values of L, used for this purpose are 0.20752, 0.797065, 1.792481 and 2.248243 for graph (a), 
(b), (c) and (d) respectively. 



V. The Middle Solution 

At this stage, we have obtained the solutions at both sides 
outside the region of the forcing terms discussed in Section 3 
and Section 4, respectively. What is left now is to 


find the solutions in the middle. This means that we need to 
solve equations (5) and (6) in the region of the forcing terms. 
We will show that the solution in this region can expressed in 
terms of Weierstrass elliptic functions. The solutions we 
obtained in Section 3, 4, and 5 will be combined to give the 
full solutions. This will be presented in Section 7. 


In this region, 


a 

2 ’ 


the solution must satisfy 


Xrj + arf + = 1 , X > 0 (23) 
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We shall be dealing with equation (23). In order to solve this equation we need to deduce some conditions. 

i a a 

These conditions come from continuity of Tj and Tj at X = — — and X = —. 


•HHO 


G and 

4i) = ' 7 ' 

(24) 

- — seek 2 
2a 

1 — — (x + L ), 

\ A/3 ^ °’ 

(25) 

T 

— rj n tanh 

P \ 


(26) 


At X = — 
2 

we have 

■ -U) 

-If) 

(26*) 


Integrating equation (23) from -a/2 to x (< a/2), we find 


Equation (27) can be written a 


=rj--r? 2 -^V 3 +Vo ( 27 ) 


[tj' ) 2 = b{q 3 + b 2 r ) 2 + b 2 rj + b 4 (28) 


By making the transformation 


b, =- —, b z =- — , b,=— andb,=b,r/ n 

1 3 p 2 J3 J3 4 3 /0 


Tj = CjM +C 2 

Equation (28) is converted into 

Cj 2 m /2 =b 1 (c l u + c 2 ) 3 +b 2 (c l u + c 2 ) 2 +b 3 (c 1 u+c 2 ) + b 4 (29) 

Dividing both sides of equation (29) by (c } ) we obtain 
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(tt 7 ) 2 =(ViV +^^-^-\ + b 2 \u 
(byC 2 3 +b 2 c 2 2 +b 3 c 2 +b 4 )— 




-2b 2 \- t J + b 3h + 


(30) 


Equation (30) is now written in the a simplified form as 


where 


u' =4 u 3 -g 2 u-g 3 (31) 


4 

c, = —, c, 

*t 


A. p - ( b 2 C 2+b 3 ) 

07 ’02 

5b l Cy 


and 



c \ 


In equation (31), g 2 is constant and g 3 is a function of L a for given A, OC, /? and a. The general 
solution of equation (31) can be expressed in term of Weierstrass elliptic function (Gong, 1994) 


u(x) =p(x + T,g 2 ,g 3 ) 


where T is a constant. 


Thus, in the region of the non-zero 
written as 


M u 

forcing term (|x| < —) the supercritical solitary wave solution can be 


tl(x) = Cyfc)(x + T,g 2 , g 3 ) + c 2 (32) 

Equations (25) and (26) are now respectively reduced to 

do =cM~ + T, g 2 ,g 3 ) + c 2 (33) 
and 

7i =c 1 p l (-^ + T,g 2 ,g 3 ) (34) 


From Gong and Shen (1993), we find the identity 


(x + T,g 2 ,g 3 ) = - 


f p / (x,g 2 ,g 3 )-p / (r,g 2 ,g 3 ) ]" 

1 P(x,g 2 ,g 3 )-p(T,g 2 ,g 3 ) ] 


-p(x,g 2 ,g 3 )-p{T,g 2 ,g 3 ) 


(35) 


Substituting equations (33) and (34) into equation (35), we then obtain 


(x + T,g 2 ,g 3 ) = - 


P* (x,g 2 ,g 3 )~ — 

C \ 

P(x,g 2 ,g 3 )-( ? 0 ° 2 ) 

C 1 


-p(x,g 2 ,g 3 )- 17 ° Cl 

C) 


(36) 
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We now observe that in equation (36), the right hand side is independent of constant T. We then have 


(x + T,g 2 ,g 3 ) 


1 f C 1 p'(x,g 2 ,g 3 )-Th V 

4\p(x,g 2 ,g 3 )-7j o +c 2 f 


-pUgi’gs)- 770 Cl 

Cl 


(37) 


Using Mathematica program Substituting equation(37) into equation (32), we get the solution on the 
middle side. All graphs of these solutions are shown in Figure 3. 

Figure 3 shows the solutions of the sfKdV in the middle, when A = 3 and a = 2. The values of L a are 
-0.20752, -0.792481, -1.792481 and -2.248243 for graph (a), (b), (c) and (d) respectively. 



Figure 3: The Middle Solution 


VI. The Phase Shift L o 


In this section we shall determine the phase shift L a which depends on the forcing terms and the continuity 
conditions between each adjacent part. 

Substituting x = a/2 and x = -a/2 into equation (32) we get 

7 (f) = <ri# {f + r,g2,g3 ) + Ca ^ + T ^2+g^j + C l- 

Substituting the above expression in equation (26*) we obtain 

<»[|+r,g 2 ,g,j-^-|+r,* J +g,j = o (38) 

From equations (33), (34), (37) and (38) we can obtain 
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B(A,L o )^ 


[ c 1 p / (a,g 2 ,g 3 )-7? l y 
[c 1 p(a,g 2 ,g 3 )-7/ 0 +c 2 J 


-c 1 p(a,g 2 ,g 3 ) + 2c 2 -2r] o =0 


(39) 


Equation (39) determines the phase shift L 0 . 

Equation (39) have a solitary wave solution only when B(A,L o )=0 (Gong,1994). 
Hence, B(ji, L a ) = 0 is the condition used to calculate the phase shift L a . By using Mathematica 
program the function B(A, L Q ) was plotted versus trial L a . The intersection point with L g -axis was 
found as the values of L a .(see Figure 4) 

This curve has four intersection point with L a -axis, i.e. B(A, L 0 ) has four zero, L 0 = —0.20752, 
-0.797065,-1.792481 and —2.248243 that means we have four supercritical solitary wave 
solutions. 



VII. Complete Solutions 


solutions of the sfKdV equation. This means that we need to 
combine the solutions in both side with the middle ones. 


We have obtained the solutions of the sfKdV equation at both We rewrite the solutions outlined in Section 3 

sides and also in the middle as discussed in Section 3, 4, and 
5, respectively. In this Section we shall give the complete 

1 (x)= -^sech 2 F^(x-L o ), (40) 

2a y 4/? 

for the left side, and 

,W= (41) 

for the right side. 
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The solution in the middle parts 


T](x) = c t ip(x + T,g 2 , g 3 ) + c 2 (42). 


By matching equations (40), (41) and (42) we get 


3A 

2a 


seek 2 



(x-L a ), 


V(x) = 


c^(x + T,g 2 ,g 3 )+c 2 - — < x < —; 

3/1 , 2 I ~X a 

- sec h - (x-L ,), — <x<oo. 

2a \ 4 P 2 


(43) 


Equation (43) is the complete solution of sfKdV equation. 

The graphs shown in Figure 5 are the complete solutions of the sKdV equation when /L = 3 and a = 2. The 
values of L o are -0.20752, -0.797065, -1.792481 and -2.248243 for graph (a), (b), (c) and (d) 
respectively 


A- a 

|/Y 

, - 2 _! ^ ! 2 3 U -- Z ■— - 

A t\ A A ^ 

Lza. 

A?vV • ■ •• -a 


-V -a AM i « » <• -* 


Figure 5: The Complete Solution 

four different solutions according to the value of the phase 
VIII. CONCLUSION shift L o . Only one solution is positive the three others 


In this paper we have studied the solutions of the stationary 
Korteweg-deVries (sfKdV) equation. In order to obtain the 
solution of the sfKdV equation, the domain of the problem 
has been divided into three parts, the left, the middle and the 
right parts. The solutions on the left and right parts were 


solutions have positive and negative values.. Further research 
can be carried out for negative forcing terms 
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